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Abstract 

We investigate the effects of all flavor blind CP-conserving unparticle operators on 
5th force experiments, stellar cooling, supernova explosions and compare the limits 
with each other and with those obtainable from collider experiments. In general, 
astrophysical bounds are considerably stronger, however they depend strongly on the 
dimension dy of the unparticle operator. While for djj = 1, 5th force experiments yield 
exceedingly strong bounds, the bounds from stellar and supernova cooling are more 
comparable for du = 2, with stellar cooling being most restrictive. Bounds on vectorial 
unparticle couplings are generally stronger than those on scalar ones. 



1 Introduction 



Recently the possible existence of a non-trivial scale-invariant sector with a non-trivial fixed 
point was proposed by Georgi [1]. These new fields, which couple weakly to Standard Model 
(SM) particles, are quite different from other extensions of the SM as they are not described 
in terms of particles but rather by "unparticles". A different, but in effect similar deviation 
from the standard model has been proposed by Van der Bij [2]. The picture is valid up to 
a certain scale, above which the picture changes. At low energies the unparticle sector is 
characterized by a scaling dimension dy, which is in general non- integer. 

In this paper we want to assess possible effects of this extension of the standard model 
in astrophysics. There are by now several studies in this direction [3,4]; in this paper we 
combine the different manifestations and give also a more detailed and complete treatment 
of the various unparticle operators. If the conformal invariance is not broken in the infrared, 
as it is assumed throughout this paper, astrophysical constraints can highly restrict the 
interactions between unparticle and SM fields. 

We consider only couplings between SM singlet unparticles and Standard Model fields 
through CP-conserving and flavor blind interactions. In Ref. [5] a list of operators composed 
of SM fields with dimensions 4 or less has been given. For our purpose we only need the 
couplings to fermions and gauge bosons. With fermions we have: 

A*// = T^T /7/J °u + i%h f^f °u + %f!PfOu + ^ Hf d^O u (1) 

ly U 1X U ly U ly U 

- ^ Hf OZ + ^ /7,7 5 / 0£ + ^ HP! O u + ^ Hf d^O u (2) 
z z z z 

+ ^hm 5 fOu+^h,i 5 f^o u . 

Here the coefficients in eq. (T5]) have been scaled to a common mass, chosen as the Z-boson 
mass Mz, so that the only unknown quantities are the dimensionless coupling constants 
ex- In the above equations, D is the covariant derivative, which introduces four-particle 
couplings involving the photon through = <9 M + ieQA^ + . . . The term with c$2 does not 
contribute to physical processes due to current conservation. 
For photons we have 

C C~~ 

r 7J - 77 f FM" o,, 77 r ^p°F f n 7J (•*) 

~ 4A du 4A du r ^ r f> a Uu W 

C TL_ f pp- u q u ° 77 e H»P° F f Ou (4) 



4M^ M ^ 4M^ M 

The term with c^p has the same structure as the effective coupling of axions. Possible 
couplings to gluons are not considered. 

In the following sections, we will consider the various standard tests of new particles and 
forces and reach our conclusions. Since unparticle phase space integration is more involved 
than for usual particles, we have added an appendix with some useful results. 
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2 Constraints from 5th force experiments 

New massless or light degrees of freedom can mediate new forces between SM particles that 
lead to an effective modification of the Newtonian law of gravity [6,7]. The most stringent 
constraints come from composition-dependent experiments, which were originally pioneered 
by Eotvos, Pekar and Fekete [8]. They make use of the fact that a new (fifth) force would 
in general act differently on different bodies of equal mass, depending on their chemical 
composition [9]. Limits for such an interaction have been derived on different length scales 
ranging from sub-meter to astronomical scales of order AU. 

For the experimental analyses, the fifth force is typically parametrized by the potential 

—r/X 

V 5 (r) = ±iGm\ 1 H l )B l B 3 - , (5) 

r 

where the coupling constant £ has been normalized to the gravitational interaction between 
two hydrogen iH 1 atoms, with G Newton's constant of gravity. Bij are the baryon numbers 
of the two test objects. 

The potential for interaction due to vector unparticle exchange can be derived from its 
propagator [10, 11] 

a- ^ / Wo|to£(OT> = 

16tt 5 / 2 T(du + 1/2) 
du (2tt) 2 *< T{du-l)T{2d u )' {) 

By taking the Fourier transform of this propagator in the low-energy limit one obtains 

Vu = Cau^L, (8) 
a u = ^Ml~ 2du ^ Y{2d u - 2) (9) 

47T 7T 
_ c 2 ,,2-2 ^(4,-1/2) 

where C = 0(1) accounts for the quark and electron density inside the nucleons and atoms 
of the test objects. For a conservative limit, we take C > 1. 

Also ca contributes in the same way as Cy, up to a prefactor, and yields 

a - (37r)V2 c 2 M 2 - M * g^^ flirt 

The scalar and pseudo-scalar interactions oc Cs2 5 cpi,cp2 do not contribute to long-range 
non-relativistic forces. However, the contribution from csi gives 

tt 1/2 2 m tmj T(d u - 1/2) 
^-(2^ Csi M 7 M « T(d u ) ' (Uj 
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Figure 1: Limits on vector unparticle interactions from Eotvos-type fifth- force experiments 
at different length scales A, for various scaling dimensions du- The dashed lines indicate the 
overall limit derived from the whole A range. 



where rriij are the masses of the electrons and nucleons between which the interaction is 
exchanged. The major contribution here comes from the nucleons with m^j pa m(iH l ) pa 

The experimental limits on an interaction of type eq. © can be applied to the unparticle 
force eq. by observing that the constraints on eq. §5$) come mainly from measurements 
at a length scale r pa A. For r ^> A, V 5 is exponentially suppressed, while for r <C A the 
experiments are less sensitive [7]. Therefore the exclusion limit at length scale A is 

« w ,um « e-^umGm 2 ^ 1 ) A M «- 2 . (12) 

This result agrees well with the power-law analysis in Ref. [12] for du = 2. 

Taking the experimental values (see Ref. [7,13] and references therein), results are shown 
for different scaling dimensions in Fig. [TJ They can be readily translated to the axial- vector 
and scalar cases. 

3 Constraints from stellar cooling 

Constraints from stellar cooling on fermion couplings. In the hot and dense envi- 
ronment of stars, light weakly interacting particles can be produced efficiently and would 
contribute to the cooling of the star. Constraints on such particles can be derived from white 
dwarfs [14-16], the ignition condition for type I supernovae [17], horizontal-branch stars with 
a helium-burning core [18-20], and red giants near helium ignition flash [16,21-23]. These 
processes have been studied extensively for axion emission, with the strongest bounds coming 
from helium-burning stars and red giants. In the following, we will focus on the evaluation 
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Figure 2: Feynman diagrams for unparticle emission in (a) Compton-like processes, (b) 
bremsstrahlung-like processes and (c) processes with unparticle-photon couplings. 



of unparticle emission from helium-burning stars, which would lead to a reduction of the 
lifetime of the horizontal-branch stars. 

Mainly two processes contribute to energy loss from horizontal-branch stars, the Compton 
process 7 + e — > e + X and bremsstrahlung involving Hydrogen and Helium nuclei as well 
as electrons, e + H+ -> e + H+ + X, e + He 2+ 



e + He 2+ + X, e + e -> e + e + X. Here 



X is the axion or unparticle. The corresponding Feynman diagrams are shown in Fig. [2] (a) 
and (b). 

The total cross-section for axion emission through the Compton process is 

2 



3m 2 . 



m P 



(13) 



in the limit u <C m e , where u is the incoming photon energy. g aee is the axion-electron 
coupling. 

For unparticle emission, the calculations are somewhat more complicated. Due to the 
phase space factor A du 6(p1 ( )9(pl ( )(pl ( ) du ^' 2 [1], the final state integration requires some care. 
The important integrals are collected in the appendix. In the limit u <C m e one finds for 
the Compton process production of unparticles 
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In the hot environment of a star photons are generated thermically, with a distribution 



n 7 (T, a;) 



7r a; 



The thermally averaged unparticle energy emission rate is then 
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with the electron density 



l + X H p 
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where p is the total density and X-r the mass fraction of hydrogen. The averaging gives 
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The emission rate for axion bremsstrahlung from electron-nucleus collisions is, in the limit 
for small incident electron velocities $ <C m e [21,24] 
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where 77/g z ctrc the electron and nucleus densities and Z is the proton number of the nucleus. 
For unparticle emission through bremsstrahlung in the non-relativistic limit one finds 
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One can see that the rate for the axial vector vanishes for djj — > 1; below we will find this 
behaviour also for other processes. This results holds however only for the leading power 
in the velocity $ in the non-relativistic limit, as one may see when expanding the correct 
expression in powers of (3. The suppression can be understood from the fact that electron- 
nucleon scattering is independent of the chirality of the particles and therefore the L — R 
coupling is suppressed. Because however unparticles for du > 1 carry a third polarization 
degree of freedom, the suppression is not total. 

The bremsstrahlung emission rates have to be averaged over a Maxwellian distribution^ 

n>e(T, Pi) = (£^) 3/ W «P("^) ( 2 9) 

so that 

/■oo 

Q%(Pi) = C(3l =► Q%(T)= dAn e (T,A)Q^f(A) (30) 

Jo 

= C2iT- l l 2 T(^f)(2T/m c y/ 2 . 

Furthermore, summing over the relevant nuclei, 

n e V Z 2 n z « n e (n u + 4n He ) « ( -?-) ■ (31) 

Bremsstrahlung in electron-electron collisions leads to very similar results as bremsstrahlung 
in electron-nucleus collisions, except for the replacement Z 2 n e n z — > An 2 in Q{j3i) or Z 2 n e n z — > 
v^2^g in Q{T), respectively [21]. Bremsstrahlung in nucleus-nucleus collisions is negligible 
since the radiation of unparticles from nuclei with mass m z is suppressed by powers of 

The impact of weakly interacting particle emission on star cooling can be evaluated with 
a numerical code for stellar evolution [22,23]. For simplicity, we give here the comparison 
of the unparticle emission rate to the axion emission constraints which have been analyzed 
earlier [18-20]. The relation between the two is summarized in Tabled] 

One needs to observe that both bremsstrahlung and Compton processes play a role in 
red giant environments. At typical horizontal-branch star densities p ~ 0.6 x 10 4 g/cm 3 and 
temperatures T w 10 8 K = 8.6 keV, the bremsstrahlung process contributes roughly 10% of 
the total axion emission rate, while the Compton process accounts for 90% of the rate [23]. 
Then the bound g aee <2x 10~ 13 [19,20] translates into the limits in Table [2J 

Constraints from stellar cooling on photon couplings. If unparticles only couple to 
photons, they would mainly contribute to star cooling through the process 7 + e — > e + U 
via t-channel photon exchange (usually called the Primakoff process), see Fig. [2] (c). In the 

1 Since the density of horizontal-branch stars is relatively low, screening and degeneracy (Pauli blocking) 
effects are negligible [14,23]. 
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Table 1 : Comparison of unparticle emission rates to axion emission rates in a stellar plasma 
of temperature T. Separately shown are the rates from the Compton process (Q c ) and the 
hremsstrahlung process (Q eZ ), as well as different values of the scaling dimension d U - 



limit uj <C m e the cross-sections for this kind of process are 
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where k is the inverse Debye-Hiickel radius, which accounts for the screening of the Coulomb 
potential of the electron in a free stellar plasma [21,25,26]. p F q are generalized hypergeo- 
metric functions. 

In previous studies, limits have been derived for the coupling of axions to photons. By 
comparing the unparticle production cross-section to the cross-section for 7 + e — > e + a, 
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these limits can be translated to corresponding limits for the unparticle couplings. Using 
T « 10 8 K = 8.6 keV, k 2 = 7.5 x 10~ 8 GeV 2 and g ail m e < 5.5 x 1(T 14 [21], we find the 
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limits in Table [3j Note that the dependence of the results on k is very mild; changing k, 2 
by an order of magnitude changes the limits in Table [3] only by up to 20%. Therefore these 
results should be reliable even without a detailed numerical code for stellar evolution. 



4 Constraints from SN 1987A 

Unparticle emission would also influence supernova cooling. This has been analyzed for 
vector unparticles in Ref. [3,4]. Here the analysis in Ref. [4] is extended to derive limits for 
other types in unparticle couplings, as in eq. §Z§. 

The observation of the length of the neutrino burst of the supernova SN 1987A puts a 
strong constraint on the allowed energy loss rate due to unknown very weakly interacting 
(un)particles [18], 

Qx < 3 x 10 33 erg cm" 3 s" 1 . (35) 

Several processes can contribute to unparticle emission from the supernova core. The dom- 
inant effect comes from neutron bremsstrahlung, n + n^n + n + U, while proton brems- 
strahlung is less important since the proton density in supernova cores is smaller than the 
neutron density. In principle bremsstrahlung processes with electrons, e + n — > e + n+U and 
e + e— > e + e + U, can be important due to collinear enhancement. However, the collinear 
phase space region is suppressed due to strong Coulomb screening effects in the dense core 
plasma, see e.g. Ref. [21]. 

Since the supernova core temperature T « 30 MeV is much smaller than the neutron 
mass, the neutron bremsstrahlung process factorizes into a "hard" nn collision process and 
"soft" unparticle radiation from one of the external neutrons. Here one can distinguish 
between the case when the bremsstrahlung coupling is insensitive to the nucleon spin (vector 
and scalar couplings of the unparticles to the quarks) [27] and when the bremsstrahlung 
emission couples to the nucleon spin (axial- vector and pseudo-scalar couplings) [28]. 

For vector and scalar unparticle-quark interactions, one finds in the non-relativistic limit 
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where /3{ is the incident neutron velocity and (Tq 11 ~ 25 x 10~ cm is the typical nn scattering 
cross section at the given energy [27]. n n ~ 3 x 10 14 gcm -3 denotes the neutron density. 
Convolution with the Maxwellian thermal distribution gives 
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Our result for <3™ V (T) has the same dimensional dependence as in Ref. [4], but we are able to 
identify an additional numerical prefactor between 0.004 and 0.0014, depending on dy. Thus 
we arrive at somewhat weaker bounds for the unparticle interactions. In addition we obtain 
bounds for scalar interaction between unparticles and Standard Model fermions. Assuming 
C > 1, the bounds in Table [2] are obtained. 

For the emission of axial-vector and pseduo-scalar unparticles, the matrix elements fac- 
torize in a similar way into the on-shell nn collision process and soft radiation from one of 
the external legs. Since the axial- vector and pseduo-scalar unparticle emission couples to the 
spins of the nucleons, one needs to take into account the spin dependence of the nn transi- 
tion, which is given by the dynamical spin structure function [28,29]. Following Ref. [28,30], 
we obtain 
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where T a is the spin fluctuation rate. Using a one-pion exchange model for the nucleon 
scattering kernel, one obtains the estimte T a w 450 MeV for the typical temperature and 
density inside the supernova core [29]. A more robust evaluation based on experimental 
nucleon scattering data [28] finds a smaller value for the spin structure function, which can 
be parametrized by using T a w 100 MeV. Taking this value and C > 1 as before, the bounds 
in Table [2] are derived. 



5 Comparison to reach of collider experiments and conclusions 

In Tables [2] and [3] we summarize our limits on unparticle couplings derived from astrophysical 
constraints. The bounds correspond to the 90% CL experimental error of the astrophysical 
observation^], for the case that only one of the unparticle couplings cx is non-zero at a time. 
For comparison we also show earlier results for limits from current (LEP, Tevatron) and 
future colliders (LHC, ILC). To get an estimate of the possible reach of a future international 
linear collider (ILC), we have assumed that it can perform the same kind of measurements 
as LEP, but with a 1000 times higher luminosity. Of course, only a proper analysis can go 
beyond this order-of-magnitude assessment of the sensitivity of ILC. The blanks in the table 

2 In the derivation of the stellar energy loss constraints, large systematic uncertainties could arise in the 
calculation of nuclear interactions and stellar evolution. Since these errors are difficult to quantify they have 
not been taken into account here. 
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indicate that no results are available from the literature for the given interaction. Some of 
the processes are not sensitive to a certain coupling, as denoted by a bar in the table. 

It can be seen that the constraints for astrophysics are generally considerably stronger 
than those from colliders. The strongest bounds are for vector/axial couplings. For small d u 
limits on a 5th force are by far the dominant constraints; however for du tending towards two 
all constraints become similarly important; here star cooling provides the strongest bound. 
For scalar and pseudoscalar couplings the bounds are generally weaker, which is mainly due 
to the higher dimensionality of the interaction operators. For du — 1, the unparticle scaling 
behavior corresponds to a regular massless particle, so that our limits also apply to any 
model which includes a new massless scalar or vector particle (see also Ref. [31]). 

For the unparticle-photon couplings, our bounds from star cooling are much stronger 
than the limits from supernova cooling, taken from Ref. [32]. These couplings could also be 
constrained by the process e + e~ — > 7 + U at LEP and ILC, but this has not been analyzed 
so far. 

This analysis is restricted to the leading CP-conserving and flavor-diagonal unparticle 
interactions. The astrophysical constraints are not sensitive to operators that involve flavor 
changing neutral currents, which can be tested in precision experiments at low energies, 
such as heavy-flavor mixing and decays [39-45], as well as to operators that only couple to 
third-generation fermions [46,47], heavy gauge bosons [48] or the Higgs boson [37,49,50]. 
Furthermore, direct CP-violation in the unparticle operators [51] can lead to new effects, 
which cannot be tested in astrophysics. 
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Energy loss from stars 


5- KT 15 


2.5 • 10~ 12 


1 • 10~ 9 


3.5- 10~ 7 


6.3- 10" 15 


2-10- 12 


7.3 • 10~ 10 


3- IO" 7 




SN 1987A 


1 • IO" 9 


3.5- 10~ 8 


1 ■ 10~ 6 


3 • IO" 5 


2- lO" 11 


5.5 -10- 10 


1.5 • 10~ 8 


4.1 • 10~ 7 




LEP 


0.005 


0.045 


0.04 


0.01 


0.1 


0.045 


0.04 


0.008 


[33, 34] 


Tevatron 




0.4 


0.05 












[35] 


ILC 


1.6 • 10~ 4 


1.4- IO" 3 


1.3 ■ 10~ 3 


3.2 • 10~ 4 


3.2 • IO" 3 


1.4 • 10~ 3 


1.3 • IO" 3 


2.5 • 10" 4 




LHC 




0.25 


0.02 












[35] 


Electroweak precision 


1 


0.2 


0.025 




1 


0.15 


0.01 




[36] 


Quarkonia 




0.01 


0.1 


0.45 










[37] 


Positronium 




0.25 








2 • 10~ 13 


2 • 10~ 8 


0.03 


[38] 



Coupling 


csi 


Cpi, 2c P2 


du 


1 


4/3 


5/3 


2 


1 


4/3 


5/3 


2 


5th force ("Eotvos") 


6.5 • IO" 22 


1.2 • 10~ 13 


1.6 • 10~ 8 


1.7- 10~ 3 










Energy loss from stars 


1.3 • 10~ 9 


7- IO" 7 


3- 10" 4 


0.13 


4- 10~ 8 


1.1 • IO" 5 


3.3 • IO" 3 


1 


SN 1987A 


8- 10~ 8 


2.4- 10~ 6 


6.6- 10~ 5 


2- 10~ 3 


5.5- IO" 8 


1.3 • IO" 6 


3.5 • 10~ 5 


9- IO" 4 


LEP 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


ILC 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 


> 1 



Table 2: Comparison of limits for unparticle-fermion couplings from astrophysical constraints and from present and future 
collider experiments. The astrophysical bounds have been derived in this work, while the collider bounds have been taken 
from the literature, as indicated by the references in the right column. Blank spaces are left where no results are available 
from the literature, while the bars denote that no bound on the coupling can be determined. 



Coupling 




du 


1 


4/3 


5/3 


2 


Energy loss from stars 
SN 1987A 


5.5- 1(T 14 
9- KT 7 


1.7- 1CT 11 
4- KT 6 


5.3 • 1(T 9 
4- 1(T 5 


1.7- 1(T 6 

8 • lir 4 



[32] 



Table 3: Comparison of limits for unparticle-photon couplings from astrophysical con- 
straints. The bounds from star cooling have been derived in this work, while the supernova 
bounds have been taken from the literature, as indicated by the reference in the right column. 



Note added 

Shortly before finishing this manuscript, we became aware of related work on 5th force 
experiments [52] where similar, though weaker limits were obtained, since these authors 
included only results from Newtonian-law experiments at short but not at astronomical 
distances. 
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Appendix: Phase-space integrals 

In the following the relevant phase-space integrals for the unparticle emission processes in 
this article are summarized. 

For the Compton process e(p) + j(k) — > e(p') +U(k') in the non-relativistic limit, with 
the initial photon energy A; = oo <C m c it is useful to choose a reference frame where the 
electron in the initial state is at rest. The phase space integration then yields 



A du f dV 1 f d 4 k' 



' - I >l 1 ' 1 Vl " 6{k' )6{k' 2 ){k' 2 ) d "- 2 {2n) 4 5 i4 \k' +p' -k-p)\M\ 2 



u 4m c u J (27r) 3 2m c> / (2vr) 4 

pi p2ui cos 9 p/ 

/ dcosiV / dp'p' 2 (2up'cos9 p , -p' 2 ) d "- 2 \M\ 2 , (43) 

Jo Jo 



12 \ c 



where \A4\ is the squared and spin-averaged matrix element, and 9 p / is the angle between 
the incident photon and the outgoing electron. The following integrals appear: 



I 



2u) cos 6i 



r(2rf^ + n) 

(44) 

f d cos 9 p , (cos P O Mw+n = 9 , + , (45) 
Jo ZcL u + n + i 
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with n — 1, 2, 3, . . . 

For bremsstrahlung e{p) + Z(q) — > e(p') + Z(q') +U{k') one finds in the non-relativistic 
limit 

au -Am e m z pJ (2^ 2m c y (2tt)3 2m J (2^^°^ H * j 

X (2tt)V 4) (A;' + p' + q' ~p-q)\M\ 2 



x(W(tf-(3 2 f r-\kT) d "- 2 \M\ 2 , (46) 



where fiij are the velocity of the incoming and outgoing electron, respectively, $ = \p\/po, 
flf = WI/p'oi and m z i s ^ ne mass °f the nucleus. After including the matrix element, the 
typical integrals are 



I 



d\k'\\kT(irn 2 M-(3j) 2 -\kf) d -- 2 

±n\ 

2_l ( r f{(3 2_ (3 2 )) 2 du+ n-^ „ = 1 >2 ,3, . . . , (47) 



r(^/- i)iy 

2T(d w + 2 



n— 1 ' 



[k'-{p-n 2 = W\ 2 {p-it)\ (48) 



/ 
/ 



1 = ^_log^±^ (49 ) 
dfip/ 1 2 



4vr (jf-j7)< m *(fi-0})*' 



(50) 



; ^ w - ^ = - 2 ; ( ^ r g^ } -i,2,3,..., (5i) 

^ dB Ft (B 2 8 2 ) 2d » In- ^l±h. - t!l csc ( 2?r ^) oU u+ 2 ( , 9] 

o WfPAPi-P,) log ^ - 2(2dM + 1) r ( -2^)r ( 2^ + 3/2) & • (52) 

Most of the above integrals are valid only for du > 1. 

For bremsstrahlung off a neutron pair, n(p) +n(g) — > n(p') +n(g') + U(k'), the situation 
is very similar to electron-nucleus-bremsstrahlung, albeit there are some differences due to 
fact that this process can only be evaluated by factorizing the strongly interacting nn — > nn 
scattering. In the center-of-mass frame 



°u _ A d u I d V / dV / d 4 A;' 6(k' )9(k' 2 ) (k' 2 ) d "- 2 5^ (k'+p' + q'-p- q) \M | 



2 



J d 3 p > J d 3 q > (2tt) 4 (5( 4 )(p / + q'-p-q) \M 



2 



2 



X 



m 2 M ~ Pj) 2 - \kf) d «- 2 \M\ 2 /\M \ 2 , (53) 
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where Ado is the spin-averaged squared matrix element for nn scattering, which in the non- 
relativistic limit does not depend on the kinematic variables of the external particles, so 
that the integration in the denominator is trivial. As before, (3ij are the velocities of the 
incoming and outgoing neutrons, respectively. 

Besides the integrals eqs. (|47|) [with m c /2 — > m n ] and (!5T|) one needs the integrals 



To arrive at the cross-section formulae in sections [3] and HI relations between T-functions 
have been used extensively in some cases. 
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